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We perform extended numerical simulation of the dynamics of dry friction, based on a model 
derived from the phenomenological description proposed by T. Baumberger et al. In the 

case of small deviation from the steady sliding motion, the model is shown to be equivalent to the 
state- and rate-dependent friction law which was first introduced by Rice and Ruina Q on the 
basis of experiments on rocks. We obtain the dynamical phase diagram that agrees well with the 
experimental results on the paper-on-paper systems. In particular, the bifurcation between stick- 
slip and steady sliding are shown to change from a direct (supercritical) Hopf type to an inverted 
(subcritical) one as the driving velocity increases, in agreement with the experiments. 

PACS numbers: 05.45. +b, 46.30.Pa, 62.20.Hg, 91.30. Px 

I. INTRODUCTION 

It is well-known that the frictional resistance is independent of the apparent area of the sliding surface and it is 
proportional to the normal load with a proportional constant fi, which is known as the friction coefficient. Traditionally, 
the friction coefficient has two distinct values: the static friction coefficient fi s , determined from the minimum force 
needed to move a slider at rest, and the dynamic friction coefficient fid, used for the friction force when a steady 
sliding motion is established. Friction is often described using the well-known Amontons-Coulomb's laws: (i) both 
fi s and fid are independent of the apparent area of the contacting surfaces and the normal load; (ii) both fi s and /id 
depend on the shear characteristics of the contacting materials |y|3|; (Hi) for most cases, fid is appreciably lower than 
fi s . This standard picture is still widely accepted nowadays; however, in many cases, significant refinements must be 
taken into account for the explanation of the observed friction. 

Our common experiences tell us that the sliding of contact bodies subjected to a steady pulling velocity sometimes 
proceeds in an alternation of periods of rest and slide rather than moving steadily (this is the reason, for example, for 
the occurrence of squeaking noises) . This unstable motion consists of periods of a stick state followed by a sudden slip 
is known as stick-slip oscillations. The stick-slip motion occurs on both large and small scales. Understanding stick- 
slip phenomenon is not only important for many engineering applications, but also for understanding the mechanism 
of earthquakes (which is a stick-slip phenomenon at geological scale @). 

Obviously, stick-slip is caused by the variation of frictional resistance during sliding. In certain range of the velocity 
V, a velocity weakening behavior (corresponding to a decreasing function fid(V) with respect to the steady relative 
sliding velocity V of the contact surfaces) is often observed in many different materials including metals Q and 
rocks 0. In this range, which generally lies in low- velocity regime, the steady motion is unstable with respect to 
perturbation; this gives rise to stick-slip phenomena. On the other hand, the static friction coefficient fi s is generally 
found to be an increasing function with respect to the contact age of contact surfaces . This can be explained by 
the fact that the plastic relaxation of the stressed contact junctions during a stick period leads to an increasing real 
contact area, hence strengthening the contact. The mechanism for the velocity-dependence of fid and the mechanism 
for the age-dependence of fi s are not un-related; both can be understood using the concept of a memory length [Q. 

The purpose of thepaper is to show, by explicit numerical simulation, that the model based on the phenomenological 
description of Ref. 0] can be used to explain qualitatively (and sometimes even quantitatively) many interesting 
features found in experiments. The model incorporates the existing understanding of the age-dependence of the 
static friction and the velocity-dependence of the dynamic friction at low velocity, but does not make a distinction 
between the dynamic and static frictions. There is no abrupt change from "static" friction to "dynamic" friction 
as in the traditional description of friction; this is important for numerical simulation used to construct the entire 
phase diagram for the dynamics of dry friction. The rest of the paper is organized as follows. We first review the 
phenomenology of dry friction, focusing on the recent experiments by T. Baumberger et al. |l|,^|. We then show that 
the phenomenological description of Ref. [Q is equivalent to the well-known Rice-Ruina friction law on rock friction 
in the case of small deviation from the constant velocity motion || . The results from simulation of the model will be 
presented in the last section together with some concluding remarks. 
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II. PHENOMENOLOGY OF DRY FRICTION 



Following the pioneering work of Rabinowicz on metals IE), many experimental studies have been performed to 
study the low velocity friction properties of various materials. It is beyond the scope of the paper to review these 
experimental work on the dynamics of dry friction. Instead, we focus on a recent extensive experimental study of 
the dry friction dynamics of a paper-on-paper block-spring system (shown in Fig. 1) by T. Baumberger et al. |jL§|. 
Their experiments verified many known properties of dry friction and gave a rather complete picture of dry friction 
in various regimes. In particular, they explored systematically the dynamical phase diagram by varying the driving 
velocity V, the slider mass M and the spring stiffness K. The features of their experiments are summarized as follows. 

• The phase diagram in control parameter space (V, K/M) consists of two regions that can be characterized 
by stick-slip and steady sliding respectively and they are separated by a bifurcation curve. The character of 
the bifurcation changes from a direct (supercritical) Hopf type in the creep-dominated regime to an inverted 
(subcritical) one in the inertial regime as V increases. 

• In the steady sliding region, the measured dynamic friction coefficient fid exhibits velocity weakening in low- 
velocity range (< 0.1mm/ s) which can be fitted as 

Vd(V) \iowV= a v - b v m(V/Vb), (1) 

where Vo is an arbitrary velocity scale; and velocity strengthening at larger velocities which can be characterized 
as 

t-Ld(V)\ hi ghV=» d + VV. (2) 

• The static friction coefficient fi s was found to increase with the contact age t st as 

Hs(t s t) = a s + b s ln(igt). (3) 

• By introducing a characteristic memory length Dq, a relation between Eq. [l] and Eq. |^, i.e. 

lowV — I^S (Do/V) (4) 

can be established. Here Do, which is given experimentally as 0.9/im 0], can be interpreted as an average sliding 
displacement needed to move to new micro-contacts. 

These results are rather general (rock- rock friction, for example, exhibits similar features). The experiments show 
that the motion of the system at low-velocity (< 0.1mm/ s) is primarily controlled by a creep process. Based on these 
results, a phenomenological model of dry friction dynamics at low-velocity regime Q] of the paper-on-paper system 
has been proposed. Both linear j2| and nonlinear || stability analysis of the model near the bifurcation give excellent 
quantitative agreements with experiments. Beside this, the transient behavior in the steady sliding region of the 
system after setting the driving velocity to zero suddenly has been studied under the aid of the model jToj ] , where 
two-stage process has been observed. However, there is no systematic theoretical study of the entire regime covering 
the crossover from the "static" friction regime to the "dynamic" friction regime and the crossover to the high-velocity 
(inertial) regime. In this paper we show numerically that most important experimental features in various regimes 
can be reproduced in a single model. 



III. A PHENOMENOLOGICAL MODEL 



We follow the phenomenological approach used in Ref. 0. Consider the experimental setup shown in Fig. 1. The 
motion of the block under the influence of the external driving force as shown in Fig. 1, is assumed to be a thermally 
activated creeping motion in a periodic pinning potential biased by the external driving force. In the low bias regime, 
where the barrier heights of the effective potential are comparable with the thermal activation energy, the velocity of 
the slider is given by 



2 



where a is the typical distance between two potential minima, and t + (t_) is the thermal time for escaping from a 
given well into its downstream (upstream) nearest neighbor. If we let the corresponding barrier height to be U+ (U-), 
then we have 

1 uj AZ7 ± 

— = ^- ex Pl 1; (6) 

where a is the thermal activation energy, a can be written as N cr RT, where N cr is the number moles of the degrees 
of freedom involved in the creep motion, ujq is the oscillation frequency about the minimum of the effective potential. 
The amplitude AUq of the periodic pinning potential is assumed to increase with the dynamical contact age variable 
<f> defined in Ref. Q , hence the barrier heights are given by 

AU± = AU (<t>) T Fexta/2, (7) 

where F ext is the external force (in the experimental setup, it is the spring force induced by the driving velocity V). 

Although Eqs. ^| and ^ are valid strictly only for a time- and position- independent external force and a time- 
independent pinning potential, they can also be used in more general cases under the assumption that the changes in 
F ext and A[/o(</>) are so slow that an adiabatic approximation is valid. Combining the above equations we can write 
the velocity of the slider as 

• / ,\ w a . f F ext a At/o. . . 

x(t)^— 2smh{ — — }exp{ }. (8) 

Z7T Zu a 

In the case that F ext a 3> a, the above equation can be approximated as 

Femt= 2Mm + 2a 2™ 
a a luqCi 

We now proceed to derive Uo(<p) by considering the case of constant velocity motion [x = V). In this case 4> = Do/V, 
and the friction force (given in Eq. |l|) is equal to the external force F ext given in Eq. | This gives rise to 

M 9 (a,-KHV/Vo)) = 2 -^ + ^H^)- 

a a uj^a 

By using the average contact age <f>, we can rewrite the above equation as 

AU o (0) = ^Mg{a v -b v ln^- - Aln^^ + {b v + A)ln4>}, (10) 

where A — 2a/aMg. Given [To (</>), we can write the friction coefficient (which is equal to the external force divided 
by the weight Mg under the quasi-stationary approximation) in term of <fi 

(i(4>, i) = a v + b v In — ^ + A In — (11) 

When the velocity is not a constant, the contact age <f> is assumed to satisfy the following equation 

<K*) = i-ff- ( 12 ) 

By using this equation, the rate-dependence of the dynamic friction and the contact time dependence of the static 
friction can be taken into account properly. Eqs.^ and ^| form the basis for the phenomenological description of dry 
friction given in Ref. Q . 

We now show that this description of friction is in fact equivalent to the state- and rate-dependent friction law 
proposed by Rice and Ruina || in the case of small deviation from constant velocity motion. To make this connection, 
we choose 8 = ln(</>Vo / Dq) [9 is the state variable in the Rice-Ruina theory), then Eqs.pl] and 12 become 

fj,(9,x) =a v + Aln-^ + (b v + A)6, (13) 
^)=cxp{-fl + hA--f . (14) 
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In Rice-Ruina theory, the dynamical equations for the friction coefficient [i and state variable in the state- and 
rate-dependent friction law || can be expressed as [[ll] 

H{0, x) = Mo + A 1 In + B'O, (15) 

!9(t) = ~{0 + ln£}. (16) 

Consider the case that there is only a small deviation from the steady sliding state, we write 

x(t) = V + Ax, 
0(t) = -\n^- + A0. 

Expending Eqs. [l4| and [16] to first order of Ax and AO give the same dynamical equation for AO, 

d{A0) _ VAO Ax 
dt ~ ~Do Aj ' 



(17) 



x) = Asign(x) sinh 1 {- exp[-^-]}, (19) 



By choosing a v — fio, A — A' and b v + A — B' , Eqs. |l3| and [lj are exactly the same. Thus these two descriptions 
are in fact equivalent in the case where the deviation from the steady sliding motion is small. This is a strong 
indication that the phenomenological theory of dry friction discussed in Ref. is not limited to the explanation 
of experimental results on the paper-on-paper system, but a rather general theory that can be used as a basis for 
studying the dynamics involving dry friction on a range of materials. For example, it will be useful for studying rock 
friction, which is important in the study of earthquake dynamics. 

For the case where a is comparable with F ex ta (this can be thought of as in the static friction regime), we have to 
use the original equation, i.e. Eq. 8, which can be rewritten as 

F ext {4>,±) = MgAsmh- l {— exp[^^]}. (18) 

ujqo, a 

The friction coefficient is then given as (with quasi-stationary approximation) 
where 

A = «„+Mn^ + J 41nM ; (20) 

here, p, is an approximation of fi we describe in Eq. |ll|. In this static friction regime, the velocity x is very small and 
it can also change its sign, thus the equation for <f> needs to be modified. It is reasonable to simply use 

^) = 1_W (21) 

There are a few advantages of using Eqsjlj], and |2l] to describe friction force: (i) there is no need to make a 
distinction between static and dynamic frictions, thus there is no stopping condition (which is used to indicate the 
crossover from the "dynamic" friction regime to the "static" friction regime) to worry about; (n) the contact age is 
a well-defined quantity with Eq. (Hi) There is no numerical singularity in the dynamical equations for the sliding 
block when the velocity goes to zero and changes sign. 

To include the inertial regime, we have to take into account the velocity strengthening described by Eq. ^; we 
include this effect by simply adding a damping term rjx in Eq. [l9| 

IV. PHASE DIAGRAM OF DRY FRICTION DYNAMICS 

We now present the numerical results from the simulation of the block-spring system. Using the friction law 
described in the previous section, we can write down the dynamical equations of the block-spring system shown in 
Fig. 1. 
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Mx(t)=K(Vt-x)-fi(<p,x)N, (22) 

where 

fi((f>,±) = Asign(x) sinh _1 {i exp[^]} + rjx, (23) 

M = a . + Mn^ + AlnM, (24 ) 

Mt) = l-^. (25) 

We compare our result for the friction coefficient /i = (j,(Dq/V, V) with the experimental results when the sliding 
motion is steady (x = V and </> = D /V). All the coefficients in the friction law are determined experimentally for 
the paper-on-paper system except for the value of m which is determined as follows. We consider the deflection point 
at V = V* from velocity weakening to velocity strengthening; this is determined using dfi(Do/V,V)/dV \v=v* = 0. 
This leads to 



and 



b v exp^ 



A J 



2F *-/^p[f] 



(26) 



where fi = a v — b v ln(V*/Vo). Substituting the experimental values of a v — 0.369, b v — 0.0f4, Vq = lfim/s, 



V* = 1mm/ s H and A = O.Olf g into Eq. [26] gives rj = U.Os/m. 

Using the values given above, the friction coefficient fj,(Do/V,V) in steady sliding region is shown in Fig. 2, which 
shows an excellent agreement with experiments g]. The logarithmic velocity weakening behavior of fi(Do/V, V) can 
be found at low velocities (< 0.1mm/ s); in this regime, the motion of the slider is creep-dominated and the effect 
of the damping term rjV is negligible. In the high- velocity regime, h(Dq/V,V) becomes an increasing function of V 
which behaves as r/V asymptotically. Hence, the experimentally observed dynamical friction coefficient /tid(V) Q can 
be described by our simple expression of /z(Z?o/U, V) without imposing any extra condition. 

We explore the dynamical phase diagram in control parameter space (V, K/M) by systematically varying V (within 
10~ 2 to 10 5 /ito/s) and K (within 10 _1 to 10 5 N/cm) with M = 1.2kg (this is the value used in the experiment). 
The phase diagram is shown in Fig. 3. Stick-slip occurs below a bifurcation curve (K/M) C (V), and steady sliding is 
found above (or to right of) it. The dots refer to the first steady sliding motions observed when increasing K or V. 
The phase diagram consists of a creep-dominated regime which can be characterized by the time scale r cr = Dq/V 
at low velocities and an inertial regime with a characteristic time Tj„ = 2n(M/K) 1 / 2 at higher velocities. The 
dashed line indicates r cr = n n . In the creep-dominated (low- velocity) regime, the bifurcation occurs at a constant 
K/M (w 1.58 x 10 3 N/cm ■ kg). Fig. 4(a) shows the bifurcation from stick-slip to steady sliding as K is increased 
in this regime. The bifurcation in the inertial (high-velocity) regime, on the other hand, occurs at a constant V 
(« 1.28 x 10 3 /i77i/s). The bifurcation from stick-slip to steady sliding in this regime is shown in Fig. 4(b). Except 
for the finite slope of the experimental bifurcation curve which we are not able to reproduced, our theoretical phase 
diagram agrees reasonably well with the experimental phase diagram. 

To investigate the nature of the bifurcation, we have measured the amplitude and period of the stick-slip oscillation 
by using different K and V. The results are shown in Fig. 5. The amplitude and period decrease as the stiffness 
K is increased. The amplitude approaches zero continuously as K becomes larger and larger; this suggests that the 
bifurcation from stick-slip to steady sliding by increasing K is a direct Hopf bifurcation. This agrees with previous 
theoretical and numerical analysis |^!|. On the other hand, the transition from stick-slip with a finite amplitude to 
steady sliding (with zero amplitude) when increasing the pulling velocity V is rather sharp (the transition velocity 
is around 1.28 x 10 3 /ito/s). This supports that the bifurcation encountered by increasing velocity is of the inverted 
Hopf type as was suggested previously [0J|,|l^| based on experimental data (but was not understood theoretically). 
We have shown numerically that the inverted Hopf bifurcation can be obtained within our fricition model. 

V. CONCLUSION 

We have shown that the phenomenological description of dry friction proposed by Baumberger et al. is equivalent 
to the friction law of Rice and Ruina in the case of small deviation from the constant sliding motion. We extend 
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the phcnomcnological description of Baumberger et al. and construct a simple model that can be used to simulate 
dynamics of dry friction in various regimes. We have shown that, by associating the thermally activated creeping 
motion with a damping term which is significant only when the velocity is large, the model gives rise to satisfactory 
agreement with the experimental phase diagram. Except for the finite slope of the bifurcation curve observed in 
experiments, the essential features of the bifurcation in both creep-dominated and inertial regime can be reproduced 
using our simple model. We believe that the model will also be very useful for numerical study of earthquake models 
which are often modeled using block-spring systems. 
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FIG. 1. Schematic block-spring system. The slider with mass M is driven by a pulling velocity V through a spring with 
stiffness K. The displacement of the center of mass of the slider with respect to the track is x. 

FIG. 2. Friction coefficient fj, described by Eq. ^ vs the pulling velocity V when the sliding motion is steady, /i changes 
from velocity weakening to velocity strengthening at V = V*(= lmm/ s); this gives r\ — 14.0s/m. 



FIG. 3. Dynamical phase diagram in (V,K/M) space, obtained by varying K and V with M = 1.2kg. It consists of a 
creep-dominated regime which can be characterized by the time scale r cr and an inertial regime with a characteristic time n n 
(see text). The dashed line indicates r cr = n n . The bifurcation curve which is represented by dots, changes from a horizontal 
line (K/M ~ 1.58 x 10 3 N/cm ■ kg) to a vertical line (V w 1.28 x 10 3 ^im/s) at the crossover from creep to inertial motion. 

FIG. 4. (a) Time evolution of the spring elongation, when crossing the bifurcation curve in low-velocity regime. V = 5.0fim/s, 
M — 1.2Kg, and (from upper to lower curve) K = 10 2 , 10 3 , 10 4 N/cm. The lowest curve has been shifted vertically by the 
amount —1.0 x 10 -4 for the sake of clarity, (b) The transition from stick-slip to steady sliding in high-velocity regime. 
K = 10 2 N/cm, M = 1.2Kg, and (from upper to lower curve) V = 10 2 , 10 3 , 10 4 fj,m/s. The lowest and the second lowest curves 
have been shifted vertically by the amoumts —5.0 x 10 -4 and —2.0 x 10~ 4 respectively for the sake of clarity. 

FIG. 5. (a) Amplitude of the stick-slip oscillations vs the pulling velocity V at various values of K. The amplitude approaches 
zero continuously as K increases, whereas the transition from stick-slip to steady sliding as V is increased is rather sharp. (b) 
Period of the stick-slip oscillations vs the pulling velocity V at various K. The period approaches zero continuously either as 
K or V increases. 
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